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Abstract 

We discuss the reaction rate of the muon-catalysed t-t fusion. The reaction rate 
is determined as a function of the temperature using the model of "in flight" fusion. 
We especially take into account the effect of the fluctuation of the cross section due 
to the existence of the muon. The obtained reaction rate 5. Ox 10~ 3 /^s _1 is a factor of 
1CP 3 smaller than the experimental muonic cycling rate in the solid tritium target. 



1 Introduction 

In the liquid hydrogen isotopes mixture, muons assist the fusion through the formation 
of a muonic molecule, since the size of the muonic molecule is much smaller than that of 
the ordinary molecules and the fusing nuclei tend to stay closer. This whole process takes 
place at the thermal energies where the conventional measurements of the fusion cross 
section using a charged beam cannot be performed. This mechanism of the muon catalyzed 
fusion (/iCF) might provide us an unique opportunity to investigate, in a rather direct 
way, the fusion cross section, i.e., the astrophysical S"-factor, at extremely low energy. For 
this purpose we need to know quantitatively the effect that a muon shields the Coulomb 
potential between colliding nuclei pp. And then this shielding effects of muons should be 
removed from the S'-factor data, in order to asses the bare reaction rate correctly. From 
this point of view, the t-t/j, fusion could provide us an elucidating example. It has been 
investigated experimentally in the gas, liquid j2] and solid [3] targets. Especially in the latter 
experiments the fusion neutron energy spectrum has been determined and its distribution 
suggests that the fusion is followed by a sequential decay of 5 He: 

t + t + ii-> 5 He* + n->a + n + n + Q(11.33MeV), (1) 

where 5 He* is in the 3/2~ and l/2~ resonant states. The t-t fusion with muons has not 
received much attention as a candidate of an energy source in contrast to the d-t and 
the d-d /j,CF. This is partly because of the difficulty of tritium handling. Moreover it is 
because its reaction rate is expected to be much lower than the others due to the lack of 
the resonant muonic complex formation. In fact the cycling rate obtained experimentally: 
3.3±0.7 /is -1 [3] (15 /xs _1 [H E] ) of the t-t/u, is smaller than that of the d-tfi of the order 
of 100. Another distinctive difference of the d-t fi and the d-d [i reactions from the t-t fi 
reaction is that their cycling rates has target density and temperature effects, which are 
likely caused by 3-body collisions. Put another way, if the cycling rate of the t-t \i reaction 
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does not have the target density and temperature dependences, one can verify that the 
dependences originate from the formation of the resonant muonic complex. 

On the other hand, the cross section of the reaction 3 H( 3 H,2p) 4 He has been measured 
in the triton beam energy range E lab =30-300 (keV) [5J El E] • This energy range is much 
higher than thermal energies. The astrophysical 5-factor of the reaction has been studied 
theoretically by means of DWBA [10] and the generator coordinate method [11] . 

We determine the reaction rate of the t-tfiCF by considering so-called "in flight" fu- 
sion [T2] and compare it with the experimental muon cycling rate. At thermal energies, 
where the /iCF takes place, fluctuations of the cross section might play an important role. 
We investigate the influence of the fluctuations by using a semi-classical method, the con- 
strained molecular dynamics (CoMD) approach [T3J, [l^] . The molecular dynamics contains 
all possible correlations and fluctuations due to the initial conditions (events). In the CoMD, 
the constraints restrict the phase space configuration of the muon to fulfill the Heisenberg 
uncertainty principle. The results are given as an average and a variance over ensembles of 
the quantity of interest, which is determined in the simulation. We make use of the average 
and the variance of the enhancement factor of the cross section by the muon, that have 
been obtained in simulations of the d-tfj,CF pQ, and convert them into the average and the 
variance of the effective potential shift. If so, one can use the same potential shift for the 
case of the i-i/zCF, because of the isotope independence of the screening effect. We, thus, 
determine the reaction rate of the i-t/iCF as a function of the temperature, taking into 
account the effect of the fluctuation of the cross section by the presence of the muon. 

This paper is organized as follows. In Sec. [2] we derive the variance of the effective 
potential shift of the reactions between hydrogen isotopes by the muons. The reaction rate 
as a function of the temperature for the t-tfiCF is determined in Sec. EJ We discuss possible 
origins of the discrepancies between the experimental muon cycling rate and the obtained 
reaction rate in Sec. HI In Sec. [5] we summarize the paper. 



2 Reaction cross sections in the presence of muons 

The influence of the muonic degrees of freedom to the reaction cross section can be taken into 
account as an enhancement by the screening effect [T51 [T6] . We determine the enhancement 
factor at the incident center-of-mass(c.m.) energy E as 

f=^- (2) 

where cr(E), and ao(E) are the cross sections in the presence, and in the absence, respec- 
tively, of the muon. The <j(E) fluctuates depending on the dynamics of the 3 (or N)-body 
system [TJ. This fluctuation of the cross section can be written in terms of the fluctuation 
of the enhancement factor: A/^ as 

Aa(E) = a (E)Af IM , (3) 

where we have assumed that the fluctuation of the bare cross section is negligible at low 
temperature compared to the screened one. Through the molecular dynamics simulation 
we obtain this fluctuation as a variance of the enhancement factor. At the same time we 
determine the average enhancement factor: f^. In our previous study we have simulated 



2 



this enhancement of the cross section by the muon in the case of the d-t //CF reaction pp. 
We utilize the results from this simulation for the case of t-t f^CF reaction. In the left-top 
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Figure 1: (Left) The fluctuation of the enhancement divided by the average enhancement by 
the bound muon (the left-top panel) in the d-t fx reaction. The fluctuation of the potential 
shift divided by the average potential shift (the left-bottom panel), both as functions of the 
incident cm. energy. The arrows in the figure indicate the point where total energy is zero. 
The Right panels are same with left panels but for the d-d reaction with a bound electron. 

panel of Fig. [TJ the ratio Af^/f^ is shown as a function of the incident cm. energy. In the 
high energy limit the ratio approaches zero, i.e., the / M distribution becomes a 5-function 
(A/ M = 0) and the average f ^ approaches 1; there is no effective enhancement. In the 
low energy limit the ratio Af^/f^ is much larger than 1; this fact implies that the system 
exhibits a sensitive dependence of the dynamics on the initial conditions, i.e., the muonic 
motion becomes chaotic. The energy dependence of the ratio Af^/ is approximated well 
by the function 2.05xE~^ c 52 , where E inc is in units of keV. This curve is shown by the 
dotted line in the left-top panel in Fig. [IJ 

We write down the enhancement factor in terms of a constant shift of the potential 
barrier. Here we have assumed that the enhancement is represented in terms of a constant 
shift of the potential barrier due to the presence of the muon. The average for the dt[i 
is in good agreement with the exact adiabatic limit with the screening potential: 

UjAD) = BEt _ BEsUe „ 8.3keV } (4) 

where BE t and BE^e are the binding energies of muonic tritium and muonic 5 He ion, 
respectively [lj. We, therefore, assume that the average potential shift is equivalent to 
U^ AD \ From the commonly used expression of the bare cross section: 

a (E) = ^l e -^ E \ (5) 
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where S(E) and r)(E) are the astrophysical 5-factor and Sommerfeld parameter. The cross 
section in the presence of the muon is expressed by 

<E) = ^ZW^ {E+U ^ (6) 

in terms of the potential shift U^. Taking the derivative of the potential shift: 

Act ct 



AU^ E + U^ 



{*T,{E + tg - 1) . (7) 



Hereafter we substitute the potential shift in Eq. (JTj) by its average and thus 

Act _ A/„ _ MJ„ 



O fa E + U, 



One can deduce the average potential shift {7 M , and its fluctuation AU^ from the ratio 

a./;,/./;, 

AU, _Af, E + U, 



This potential shift is independent of the difference of isotopes, so that we make use of the 
same value in both cases of the d-t and the t-t [iGF. 

In the left-bottom panel of Fig. [1] the fluctuation of the potential shift divided by the 
average potential shift is shown as a function of the incident energy. It is striking that 
the slope of the ratio AU^/U^ changes at the ionization energy of the muonic tritium, 
which we indicate by the arrow in the figure. At this incident energy the total energy of 
the system is zero. The total system is unbound at the incident energies higher than this 
point, while the 3-body system might be bound at lower energies. By contrast, as it is 
shown in the right-bottom panel in Fig. [TJ in the case of the bound electron screening the 
binding energy of the electron is much lower than the incident energy of our interest and 
the ratio AU e /U e decrease monotonically as the incident energy decreases. Again the arrow 
indicates the ionization energy of the deuterium atom. The irregular muonic motion leads 
to smaller external classical turning point [lj. As a consequence the irregularity makes the 
enhancement factor larger opposed to the result of the electron screening [T5| [T6] , where the 
irregular(chaotic) events give smaller enhancement factors. This contradiction is accounted 
for the fact that the system remains bound in the present case at low incident energies, 
while in the previous case even the lowest incident energy that was investigated is much 
higher than the binding energy of the electrons. Therefore the chaotic dynamics of the 
electrons causes to dissipate the kinetic energy between the target and the projectile and 
lowers the probability of fusion. 



3 Reaction rate 

In a liquid hydrogen tritium target at the temperature T, the velocity distribution, <f>(v), 
of a pair of colliding particles is written as the Maxwellian distribution, 

9 W e ij W 
^(E,T)dE = 4>{v,T)dv = -=— -e-^-==, (10) 

y/ir k B T \Jk B TE 
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where E and v are the relative energy and the velocity of the pair of colliding particles, in 
the present case two tritons, and is Boltzmann constant. Although the experiment in [3] 
has been performed using the solid target, in this paper we assume Eq. ( TTU1) as the velocity 
distribution for simplicity. We will reconsider the validity of this assumption afterwards. 
The reaction rate per pair of particles is given by [17] 

< av >= J a(E)v^(E,T)dE, (11) 

where cr(E) is the reaction cross section in the presence of the muon. The reaction rate at 
the liquid hydrogen density pm — 4.25 x 10 22 cm' 2 : 

A = Plh < crv > (12) 

is obtained SIS db function of the temperature. 

As we have seen in Sec. [2j the cross section in the presence of the muon is expressed 
effectively with a potential shift U^, hence we write down the reaction rate: 

\ = Plh J <y a {E + T)dE, (13) 

where <Jq{E + U^) fluctuates depending on the variance of U^. This fluctuation of the 
screening potential can be incorporated through the following equation: 

X = Plh J ME + Ur)vV{E,T)N(Ur)dEdU„ (14) 

where N{U^) is the distribution of the screening potential. As the first guess, the distribu- 
tion is likely to be a normal (Gaussian) distribution. In reality, however, we find that the 
distribution changes characteristically as a function of the incident energy. This change is 
observed at the point where the incident energy coincides with the ionization energy of the 
muon. In Fig. [2] the left panel shows a histogram of the screening potential at a relatively 
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Figure 2: The histograms of the screening potential divided by the average at a relatively 
high incident energy (the left panel) and at an incident energy lower than the ionization 
energy (the right panel). 



high incident energy E =60 (keV). The abscissa is the screening potential normalized by 
the average value £7 M =5.6 (keV). The distribution of the histogram is approximated well 
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by the Gaussian distribution, as we expected. Moreover one sees that there are some 
components with negative screening potentials. This negative component of the screening 
potential means that muons can be kicked out to an unbound state in some cases and bring 
away the relative energy of the colliding nuclei. As the incident energy is reduced, the 
distribution changes. The right panel in Fig. [2] shows one of such a situation: a histogram 
of the screening potential at the incident energy E =0.3 (keV). In this case, because the 
muon is forced to remain bound the whole process in the entrance channel, the screening 
potential cannot be negative any more. In fact, we can approximate the distribution of the 
histogram with the distribution of Gaussian Orthogonal Ensemble (GOE) rather than a 
normal distribution. We, therefore, approximate this fluctuation of the screening potential 
by Gaussian distribution around the average potential shift and with variance AU^. 



N ^=^r\--^r)- (15) 

in the energy region E > BE t and by GOE: 

in the energy region E < BE t . Using these distributions of the screening potential, we assess 
the Eq. (jllj) . We substitute the bare cross section by Eq. (jSJ) and we use the polynomial 
expression [TO] 

S(E) = 0.20 - 0.32E + 0.476E 2 (Me Vb), (17) 
where the energy E in MeV, as the astrophysical ^-factor in Eq. (jHJ). 



4 Results and discussions 

The obtained reaction rate by Eq. (I14p is shown in Fig. [3] by circles linked with a solid 
line as a function of the temperature. The circles linked with the dotted line are obtained 
by taking into account the enhancement of the cross section using the average potential 
shift Up alone, i.e., Eq. ( fl3l) substituted by U^. The triangles show the bare reaction 
rate. The circles deviate from Eq. ( fl~31) at the temperature lower than 10 8 K and do not 
show much temperature dependence from 0.1 K to 10 5 K. The deviation of the reaction 
rate by Eq. (I14p from the one by Eq. (fTSj) in the low temperature region indicates that 
the fluctuation of the enhancement factor has a strong influence on the reaction rate at 
the thermal energy. Including the fluctuation of the cross section, the reaction rate for the 
t-tfiCF reaches at 5.0 x 10~ 3 /is _1 . This is a factor of 10~ 3 smaller than the experimental 
muonic cycling rate in the solid tritium target, 3.3±0.7 /xs _1 [3], which is marked by the 
square in Fig. [3J We remark several possible explanations of this underestimation of the 
rate in the following. The first one is that the assumption of the velocity distribution of 
colliding nuclei to be the Maxwellian could be inappropriate. We should make use of the 
velocity distribution in the solid target. Otherwise the second possible explanation is that 
assuming that the procedure of the estimation of the reaction rate is correct, the result 
could mean that the actual bare ^-factor must be greater than that given by Eq. (fTTI) in 
the low energy region where the measurement has not been performed. 
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Figure 3: Reaction rate of the tt\x system at the liquid hydrogen density as a function of the 
target temperature. The circles linked with the dotted line are obtained taking into account 
the average enhancement of the cross section by the muon, alone. The circles linked with 
the solid line are obtained taking into account the fluctuation of the enhancement of the 
cross section by the muon. The triangles show the bare reaction rate. The square shows 
the experimental data of the muon cycling rate. 

5 Conclusions and Future Perspectives 

In this paper we have estimated the reaction rate in the muon-catalysed t-t fusion as a 
function of the temperature in the framework of the "in flight" fusion. We have used the 
CoMD simulations to estimate the enhancement effect of the cross section due to the muon. 
We found that the estimated reaction rate does not show temperature dependence in the 
low temperature region in contrast to the experimental data of the d-t and the d-dfiCF. 
The obtained reaction rate in the low temperature: 5.0 x 10 _3 /is _1 underestimates the 
experimental muonic cycling rate in the solid tritium target, 3.3±0.7 factor of 10 . 

This is either because our assumption that the velocity distribution is Maxwellian is not 
correct in the case with the solid target or because the S-factor, which we used to estimate 
the rate, in the low energy region is lower than the actual value. 
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